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We demonstrate the emergence of the magnetic moment and spin-vorticity coupling of chiral 
fermions in 4-dimensional Wigner functions. In linear response theory with space-time varying 
electromagnetic fields, the parity-odd part of the electric conductivity can also be derived which 
reproduces results of the one-loop and the hard-thermal or hard-dense loop. All these properties 
show that the 4-dimensional Wigner functions capture comprehensive aspects of physics for chiral 
fermions in electromagnetic fields. 


1. Introduction. Significant progresses have been made in understanding the dynamics of chiral fermions in elec¬ 
tromagnetic fields. This is particularly interesting in high energy heavy ion collisions where very strong magnetic 
fields can be created. The magnetic fields are so strong that quarks can be polarized and their momentum directions 
are parallel or anti-parallel to the magnetic field depending on quark chiralities and charges. Quarks with the same 
charge tend to move in the same direction. Any imbalance in the number of right-handed and left-handed quarks 
as a consequence of topological configurations of gauge fields may lead to such a charge separation effect which can 
be tested in experiments [l|. This is termed as the Chiral Magnetic Effect (CME) [2, yb The Chiral Vortical Effect 
(CVE) is also an accompanying effect due to rotation in a relativistic and charged fluid [1,0. The interplay of chiral 
magnetic and chiral separation effects maWead to a phenomenon called the Chiral Magnetic Wave 0 , whose vortical 
counter part is the Chiral Vortical Wave [3| ■ 

Kinetic theory is an important tool to describe these phenomena in phase space of chiral fermions. The Abelian 
Berry potential takes an important role in 3-dimensions (3D) kinetic approach in accommodation of axial anomaly 
|§-[l(]|. It has been shown that the CME, CVE and Covariant Chiral Kinetic Equation (CCKE) can be derived in 
quantum kinetic theory from the Wigner function in 4-dimensions (4D) in external electromagnetic fields (Hill. 
The 3D chiral kinetic equation [M3 can be obtained from the CCKE by integrating out the zero component of the 
4-momentum. 

In the 3D chiral kinetic equation, it has been shown that the fermion energy is shifted by the interaction energy 
of magnetic moment with the magentic field [l3| . The magnetic moment and spin of fermions have relativistic origin 
|14l - ll6l |. It is a natural conjecture that the magnetic moment should also emerge in the covariant quantum kinetic 
approach in 4D Wigner functions. In this paper, we will demonstrate the emergence of the magnetic moment as well 
as spin-vorticity coupling in the framework of covariant quantum kinetic theory based on 4D Wigner functions. We 
will also show that the parity-odd part of electric conductivity (chiral magnetic conductivity) can also be derived from 
4D Wigner functions in linear response theor y w ith space-time varying electromagnetic fields. The result reproduces 
the chiral magnetic conductivity of one loop [l7l | and hard-thermal or hard-dense loop (HTL or HDL) [TsI . [l^ under 
proper approximations [l^ . 

We adopt the same sign conventions for the fermion charge Q and the axial vector component of the Wigner function 
as in Refs. 

2. Wigner functions for chiral fermions. In a background electromagnetic field, the quantum mechanical anologue 

of a classical phase-space distribution for fermions is the gauge invariant Wigner function Wap{x,p) which satisfies 
the equation of motion (2l|, [23, — m)W{x,p) = 0, where x = (xo,x) and p = (poiP) are space-time and 

energy-momentum 4-vectors. For the constant field strength the operator is defined by = p^ -I- i-^V^ 
with = dif — QF^'^dP. The Wigner function can be decomposed in 16 independent generators of Clifford algebra, 
whose coefficients and are the scalar, pseudo-scalar, vector, axial-vector and tensor components 

of the Wigner function respectively. For massless or chiral fermions, the equations for and are decoupled from 
other components of the Wigner function. 

3. Formal solution to quantum kinetic equations. For chiral (massless) fermions, we can rewrite the equations for 
and into those for right-handed (s = -I-) and left-handed (s = —) vectors ^^{x,p), 


P^/n^x,p) = 0 , 

= 0 , 


( 1 ) 
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where are given by 

= ]^[r^{x,p)+sj^^{x,p)\. ( 2 ) 

One can derive a formal solution of in Eq. (IB b y a perturbation method in powers of and F^y. The solutions 
at the zeroth and first order were given by Ref, (ill. [T^ and can be cast into the following form, 

/(‘q)s^X,p) = pPfsSip"^), 

c/(i)si^^p) = --^QF^^PxfsSip^), (3) 

where po = u • p with being the fluid velocity, and fs are distribution functions of chiral fermions defined by 

2 

fs{x,p) = J^^[Q{po)fF{po- Ps) + Q{-po)fF{-po +Ps)]- ( 4 ) 

Here fF{y) = l/[exp(/ 32 /) + 1] is the Fermi-Dirac distribution function, and (3 = 1/T and ps are the temperature 
inverse and the chemical potential of chiral fermions with chirality s respectively. We can decompose ps into the 
scalar and pseudo-scalar parts, = p + sp^. We have used following formula, 

and = \{dvUa — d^u^), where and e^j/o -/3 are anti-symmetric tensors with = 1 (— 1 ) and e^vrrp = — 1 ( 1 ) 

for even (odd) permutations of indices 0123, so we have _ —60123 = 1- Instead of Cxi/ and F^^, 

we will also use the vorticity vector = ^eP'^°‘^u^daUp, the electric field = F^'^u^, and the magnetic field 
= ^€>^''^Pu^Fxp. 

4 . Magnetic moment and energy shift. In this section we will derive the magnetic moment and energy shift from 
the formal solution of the Wigner function (|3]). To this end we assume that the chemical potential ps depends on 
space-time and 3-momentum. For convenience, we will focus on the electromagnetic term in this section and turn off 
the vortical term. The vortical term will be investigated separately in the next section. 

If Ps does not depend on the 3-momentum, from Eq. ©, we obtain the fermion number current jP (the electric 
current should be QjP), 

fsix) = J d‘^pfP{x,p)=nf'^uP + a^psB^, ( 5 ) 


where ax}s = is the parity-odd part of the conductivity in the static limit for chiral fermions with chirality s. 

Since fs{x,p) does not depend on 3-momentum, the momentum integral involving is non-vanishing only along 

the fluid velocity uP , so one can just make replacement pP —>■ pquP. Here nP is the fermion number density of a 
non-interacting gas of chiral fermions. For simplicity of notation, we will work in the local rest frame of a fluid cell 
with uP = (1,0,0,0), then we have Ep = |p|. We will use the velocity on-shell vector yP = (l,v) with the 3-velocity 
V = p/Ep. 

In order to extract the magnetic moment and its energy shift, we can extend our previous scenario by assuming that 
Ps depends on the 3-momentum since the energy shift from the magnetic moment can be grouped into the effective 
chemical potential, while the magnetic moment of a chiral fermion is proportional to its spin which is polarized 
along its momentum. We can write the chemical potential into the sum of a constant part and a (x, p)-dependent 
part, /r®(x, p) fv ps + Wf{x,p)^ where ps do not depend on x and p, and e = ±1 denote positive/negative energy 
corresponding to fermions (e = -I-) and anti-fermions (e = —). The quantities Wf{x,p) can be further decomposed as 
Wf{x,p) = W(x,v) - 1 - ^^W 5 {x,v), where the introduction of the factor l/(2£'p) in the second term is to make W 5 
depends on v only. Note that our definitions of /rs, W and W 5 can be matched to those in Ref. up to constants. 
In this Letter, we will not use the specific form of ^s(x,p) except its formal dependence on e, s, x and p. 

With ^®(x,p), fs{x,p) in Eq. ([4]) now become 


fs{x,p) 


2 


^ 0(epo)/F[epo - ep^x, ep)]. 

e^±l 


( 6 ) 


Note that there is a minus sign (e = — 1) in front of p in the anti-fermion sector. So the fermion number 4-current 
jP(x) can be obtained by integrating ^f{x,p) in Eq. ([3]) over 4-momentum with fs given by Eq. ([3]). The fermion 
number density can be derived in Appendix and be put into a compact form 


Us 




df'p 1 


sQ(v-B) ^ efF[Ep 

e^±l 


(^pI{x,p)\ 
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'P e=±l 


/ ^ e/F[£;p-eAi®(x,p)], 


e=±l 


(7) 


where = (l + sQa^-B) is the phase-space measure with the Berry curvature ~ Ep(l — sQa^-B) 

is the effective energy. So from Ep we can read out the energy shift and magnetic moment of chiral fermions, 


AEb 

Mm 


-Mm ■ B 



sQ 


( 8 ) 

(9) 


where S = ^esv and 5 = 2 are the spin and g-factor of charged chiral fermions respectively. For e = s = -I- or 
e = s = — (positive energy and right-hand chirality or negative energy and left-hand chirality), S = (positive 
helicity), otherwise S = — (negative helicity). However only depends on chirality and charge of positive energy 
fermions. For s = Q — + or s = Q = — (positive chirality and positive charge or negative chirality and negative 
charge), is parallel to p, otherwise it is anti-parallel to p. We also see the emergence of the phase-space measure 
^^ 7 , which is not surprising since we have already shown this in Ref. [a. What is new is the emergence of the 
magnetic moment and energy shift in our Lorentz covariant kinetic approach although superficially there is no such 
terms in the distribution function fs in Eq. We then reproduce the results of Son and Yamamoto (T^ . [^ . . 

5. Energy shift from spin-vorticity coupling. We now consider the vorticity term in the Wigner function in 

Eq. ®. The fermion number current from vorticity term is given by 



where the first term (^ ujP) gives the CVE coefficient, while the second term ^ provides the correction to the 
fermion number density from the spin-vorticity coupling energy. Note that these two terms are orthogonal to each 
other since u • w = 0, so there is no correction to the CVE coefficient at this level from the second term. We note that 
there are a few arguments about the temperature dependent part of the CVE coefficient such as its origin and whether 
there are corrections from high order terms (23 - [^ . The part of in Eq. (fTUl) gives the vorticity contribution 
to Us in Eq. ©■ So the spin-vorticity energy shift can be read out from Eq. m, 

=-ies(a; • v) =-w • S, (11) 


where uj is the vorticity 3-vector. 




(l)s 


From the energy density one can also derive the correction from the vorticity term 


£*(1) 


/ 


d^p 


e=±l 




eh't{x,p)]- 


( 12 ) 


One can also read out the same energy shift AE^^ from the spin-vorticity coupling as (ED. 

Combining Eq. ([5D and (ITTp , we obtain the total effective energy with energy shifts from the magnetic moment and 
the spin-vorticity coupling, 


E; = Ep - • B - u; • S. (13) 

We see the emergence of the magnetic moment and the spin-vorticity coupling from the Wigner function solution ([SD . 

6. Parity-odd electric conductivity. We can calculate the parity-odd part of the electric conductivity or chiral 
magnetic conductivity, k), from the Wigner function solution ([3D as the result of expansion in 

where we assume that d!f - (w , k) is of the same order as This is valid for small uj and |k|. For finite values of 
oj and |k|, a more rigorous way is to carry out the expansion of linear response in space-time varying fields only and 
regard d!f ^ (w, k) as 0(1) quantities. In this case, one should solve a different set of equations with similar structure 
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to Eq. Q but with replacements p>^ ^ pf^ - iQji (iA) Ff^^dP and ^ - Qjo (iA) F^'^dl, where A = ■ dp 
and dx acts only on F^^ but not on distributions. Here jo{y) = sm(y)/y and ji{y) = {siny — y cosy)/y^ are spherical 
Bessel functions. The solution to the leading equations is still fs{p)S{p'^), same as in Eqs. pi4l) except 

that fs{p) does not depend on x because should satisfy = 0. We refer the readers to Appendix [Bl for 

the details of Wigner equations and its solutions in this case. 

As shown in Appendix [B] the parity-odd part of the Wigner functions in fc = (w, k) representation which linearly 
depends on fields can be written in the following form 




. sQ 

*2p • p' ^ 2^ I 


. sQ 
2p ■ k 


t,,uoakPp''AP{k)jo ( -A ) (fc • dp)[fsS{p^)] 

2" 




k''pPAUfAp+lk]S 


P+lk 


-fs{p - ^k] 6 


p-^k 


(14) 


where we have used A = —ik ■ dp, jo A) {k ■ dp) = 2i sin A) and exp (^fc • dp) fsSijP') = fs{p + \k) S 
We can obtain the 3-vector current by integration over 4-momentum for the spatial component of 
— J d^p ^( 1 ) We can then read out k) which is just the one-loop result, Eq. (36) of Ref. |17l |. 


e- js(w,k) = 


cr^(w,k) 




(2|p| - wjln 


(^-IpI)^ 


+ (2|p| + Lu) In 


(C^-|P|)2- 

(^+|P|)^ 


(^+|P|)^ 


(|p| + |k|)^ 
(|p|-|k|)2 
-(|p| + |k|)^ ' 
-(|p|-|k|)^ ■ 


(15) 


If we assume that external frequency and momentum are much smaller than internal momentum, w, |k| <C |p|, we can 
reproduce the HTL or HDL result [l3, , 


HTL/HDL(^^k) = 4°) 1- 


|k|^ 


1 _ _^i ^+1*^1 

2 |k| ’^w-|k| 


(16) 


where we have used = ^^QP 5 - By including finite damping rate into fermionic propagators the 

non-commutativity of the static limit uj —>■ 0 and |k| —>■ 0 can be removed (^ . 

7. Summary. We have investigated properties of chiral fermions in electromagnetic fields in covariant quantum 
kinetic theory based on 4D Wigner functions. We have shown that the energy shifts of chiral fermions from the 
magnetic moment and the spin-vorticity coupling emerge from the Wigner function solutions derived from our previous 
works. We have also calculated the parity-odd electric conductivity in linear response theory and reproduced the results 
of one-loop and HTL/HDL. All these properties, together with previous findings, show that the 4D Wigner functions 
capture comprehensive aspects of physics for chiral fermions in electromagnetic fields. 
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Appendix A: Derivation of Eq. ([^I 

In the case with momentum dependent chemical potential p'^{x,p), the zeroth order term of which corresponds 
to Bq. @ can be derived as 

d(o)s(^) = j d‘^PMx,p)pP6{p'^) 
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/ 


d?p 

(Pp 


[{uP + vP)fF{Ep - pt (x, p)) + (-w^ + V^)fF{Ep + p, {x, -p))] 
v'' [fF{Ep - p+ {x,Tp)) - fF{Ep + p-[x.p))] . 


(Al) 


In the above equation, we have substituted fs{x,p) of Eq. ([S]), used (5(p^) = 2 ^['^(Po — Ep) + 6{po + Ep)], and carried 
out the integral over pq. We have used the notation = (0, v) with v = p/Ep = p/|p|- From the second to the third 
equality, we have changed the integral variable p —>■ — p for the second term (the negative energy sector) inside the 
square brackets. 

The first order term of jP corresponding to the electromagnetic part of hr Eq. ([3]) can be written as 


= sQpP^ J (Ppfs{x,p)px5'{p^) 

= f <Ppfs{x,p)^^5{p^) 

2 J Po apo 



d p— 
dpo 

fAx,p)— 

1 Po] 

5{p^) 


d p— 
dpo 

fs{x,p) — 6(p^) 

Po 


f dS 

PP 

dfs{x,p) 

- fs{x 

2 j 


[po 

dpo 


(A2) 


In the first equality we have used the identities S'ijP) = -^5{p^) = and = ^^(P^)- From the 

third to the fourth equality, we have dropped the integral of complete derivative since fs{x,p)^6{p^) —^ 0 at the 
boundaries po —>■ ±c». In the last equality we have used the notation p^ = (0,p). In order to further simplify Eq. 
(IA2I) . we will use 


dfs{x,p) 

dpo 


(27r)3 


0(po):^/f(po - P^(a;,p)) 
dpo 


-Q(-po) 


di-po) 


fpi-Po+Ps (a;,-p)) 


(A3) 


where we have neglected derivatives of 0(±po) because they vanish when combining with S{p^). Inserting Eqs. (I6IA3I) 
into the last equality of Eq. (jA2p . we obtain 


+£QFP>y 


d^p 


{2fYEp 

d^p 


_ djF{Ep- p+{x,p)) , , , , ,d/F(Ep + p, (a:,-p)) 

{u\ + V\) -—-h {-U\ + V\) - 

fi: • X fr) 


dEn 


V\ [friEp - pt{x,p)) + fF{Ep + p, {x, -p))] 


= - 2 O 


/ 


d^p 


(27r)3Ep 

d^p 


pp^ 


v\ 


dfFjEp - p+(x,p)) _ dfFjEp+p^ (a:,p)) 
dEp dEp 

PP^vx [fpiEp - pt{x,p)) - fF{Ep + pj{x,p))] , 


(A4) 


where we have changed in the second equality the integral variable p —>■ — p for the second terms (sectors of negative 
energy) inside square brackets, and we have also used v\ = u\ + v\. 

From Eqs. (|A1IA4I) . we can obtain the fermion number density by contraction of the fluid velocity Up and jf = 
d{o)s + = XLpjP{x), which gives Eq. d?]) by using UpvP = 0 and UppP^vx = UppP^vx = B • v. 


Appendix B: Wigner functions in space-time varying electromagnetic fields 


We have discussed in section 2 and 3 that Eq. (HD holds for the constant field strength For space-time varying 
Fpy, it takes the general form (see Eqs. (5.12-5.21) of Ref. 0), 

= 0, 
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= 0, 

2s{W-li'' (Bl) 
The operators 11^ and now replace and in Eq. o respectively, 

T\P = pP- iji Qa) QFP'^dl, 

GP = dP-3^{^l^QFP''dl, (B2) 

where A = dx ■ dp, and jo{x) = sinx/x and ji(x) = (sinx — xcosx)/x^ are spherical Bessel functions. Note that dx 
acts only on F^'^ but not on other functions to its right. For constant we recover 11^ = p^ and G^ = by 
using jo(0) = 1 and ji(0) = 0. 

The solution to Eq. (EH) can be found by perturbation in fields. To the first order, we can formally write 



The zeroth order equations read 


= 0, 

25r^(o).-PVro)J = (B4) 

whose solution is in the same form as the first line of Eq. (|31) except that fg does not depend on x here due to 
dP ~ order equations read. 


- JO QFP’'dPJ^,^^ 

- JO (^a) QFP>'di/^,^g] 


0 , 


0 , 


Pva^{l)p\ 


+s 


Jl ( 2 ^] * 3 -Pl'cr 3 p ^( 0 )^ 


Jl ( 2^ ) c/'{0)u 


Contracting d^ with the third equation and using the second equation of Eq. 

P-9x/(i)p{x,p) = ppjo[^l^QFP’^dl/(,^)p 


, we arrive at 


_ B’' 

2 


jo{^A)QFP’'dl^l,^g 


-\Qd: 


Jl f 2 ( 0 )m ^pcrdp ^(o)j 


(B5) 


(B6) 


We can solve in momentum space by replacing dx —>■ —ik and A — ?> -ik-dp. After a lengthy but straightforward 

calculation, we obtain 


f{i)pik,p) = • k){A ■ dp) - {p ■ A){k ■ dp)]jo (^^A^ [/s<5(p^)] 

(^a) (fc . dp)[fg6{p^)] 

+ - k^Apjik ■ dp)jo Qa) [fgSip^]] 

+i^Q[kp{A ■ dp) - Ap{k ■ dp)]ji Qa^ [fsSip'^)]. 


(B 7 ) 
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The second term is the parity-odd part and can be put into the form of Eq. ca¬ 
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